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ABSTRACT 

Analyses of electric power systems ordinarily include only a 
small portion of the entire network. Effects external to the network 
of interest are modeled as loads or generators. A m.ethod is proposed 
for identifyin,g an equivalent model of the external network without 
taking any measurements in the external network itself. The method is 
simple and could easily be implem^ented to supplement state estim.ation 
techniques. 

A linear load flow model which relates real power to voltage phase 
angle is developed for analyzing the entire network and from this two 
models are developed for relating m.easurements within the system of 
interest to the parameters of the equivalent m.odel. It is found that 
there are two components to the equivalent m.odel. One is a fixed struc- 
ture and the other is an equivalent power. Since both models for 
identifying the equivalent network are in the fonn of an input, output 
relation with an unknown additive disturbance, a general solution for 
identifying linear static systemjs is found subject to certain conditions 
Applying the solution to the two models for identifying the equivalent 
network, it is found that one of the m.odels has an inherent tendency to 
produce consistently biased estimates of the equivalent model paramieters 
However, that sam.e model has the advantage that it uses data for its 
input, output that is readily available and relatively accurate. 

A sim.ulation made to test the proposed method is discussed for a 
sample system and the results agree well in both the actual param.eters 
identified and predicted error. Several important points are discussed 
with a view toward implem.entation. 
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1. irTRGDUCTION 

Analyses of electrical power transmission networks usually carry 
the implication that the network under consideration can be isolated 
from all other networks. In actuality this can seldom be realized 
since the network is interconnected with one or more outside networks, 
and whatever happens to the network being analyzed' is determined at 
least to some extent on how it interacts with the others. Therefore, 
it would be convenient if the effects of all the outside networks on 
the network of interest could be summarized by one equivalent network. 
In most cases the equivalent network would be considerably smaller than 
all the actual outside networks that it would represent. For instance, 
if the Ifew England network were interconnected with the rest of the 
United States at a half dozen points, as far as the New England network 
would be concerned, the entire rest of the United States could be 
shrunk to an equivalent network having a half dozen buses. A method is 
proposed for findir^ the equivalent network which could be used for 
contingency planning and state estimation techniques C3 3, t 4 ts 3j 
rs],r93, and [ 123 . 
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1.1 General Problem; Identification of Linear. Static Svsteins 

When analyzing a systeE, the visual approach is to surround the 
system of interest with a conceptual boundary. That within the boundaiy 
is idealized in some respect to emphasize certain features. Further, 
the idealized system is allowed to communicate with the outside only 
through certain conceptual ports. This idealized representation of 
actuality is called the system model and by analyzing the model conclu- 
sioriS can be drawn about the actual system. If the contents of the 
idealized system and what enters through the ports are completely known, 
then what happens to the idealized system can be found at least in 
principle. If the model is static, then only what currently enters is 
necessary to know the current state of the system. In such a case, if 
only knowledge of the current state is required, then the effect of the 
system outside the boundary is summarized by what enters throvigh the 
ports. However, to know what will happen in the future, it is necessary 
to know what will enter in the future, and it may be that what enters 
through the ports is at least partly determined by how the system inter- 
acts with the outside. If everything outside the system is also ideal- 
ized, then everything of interest can be modeled by the two idealized sys- 
tems. The first, the original system model, represents what is within 
the boundary. The second, the external system model, represents all 
that is not represented by the first. Further, the two models can only 
interact through the originally specified ports. The external system 
model m.ay be rather extensive and it might be desirable to find an 
equivalent external model which would interact with the original system 
model in the sam;e way. 
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The process of firdirig the equivalent system model, called identi- 
fication here, is investigated for a class of linear, static system 
models. Kore specifically, it is assumed that the model consists of a 
structure which remains constant and two types of variables, inputs (or 
causes), and outputs (or effects). The mathematical notation using 
matrices is simply 2 = C u where u represents an input vector, y repre- 
sents an output vector, and C represents the structure. 

If the original system model with its conceptual ports to the 
outside is placed in this framework, then it might be possible to 
separate what exists at or passes through 1he ports into the same two 
types of variables. For the model of a physical system, the product of 
the two variables frequently represents a generalized power passing 
through the ports. An example would be the voltage level existing at 
an electrical terminal and the current passing through the tenninal. 

By combining inputs from the ports with Inputs internal to the system., 
the outputs can be evaluated. In the same way knowledge of inputs at 
the ports, inputs within the external system model, and the structure of 
the external model allow outputs of the external model to be found. That 
is, each model may be analyzed independently of the other by knowing 
what passes through the ports. However, in order to predict future 
outputs of the original system, knowledge of how the two m.odels interact 
is necessary. That is the purpose of identifying the equivalent external 
system model. 

The inputs and outputs can be divided into those belonging to the 
original system, those of the external system, and those common to both. 
Two general approaches to identifying the equivalent model are investi- 
gated. Orje involves combining both the original and external models 
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and then identifying the new model using all the inputs not in the 
external model and only the outputs common to both. Another approach 
is to use the inputs and outputs at the ports or terminals to identify 
the equivalent external system separately. Both approaches are analo- 
gous to a black box which contains both the external model structure 
and inputs belonging to the external model. Inputs from the original 
system enter the black box and outputs are those which are commion to 



both. 
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1.2 Specific Problem; Electrical Power Transmission fetworks 

Ordinarily when analyzing electrical power transmission systems, an 
arbitrary boundary is drawn about some part of the entire network and 
that part of the system is studied assuming power flowing in through 
the transmission lines on the boundary is known. The specific problem 
investigated is that of determining how an interconnected electrical 
power system affects a particular subsystem. The subsystem will be 
named "Own System” and referred to as CS for convenience. The rest will 
be named "External System.” and referred to as XS. Together CS and XS 
form the entire system named ”T,i/hole System” or WS . The External System 
may be the entire outside world as viewed from one region. Or CS may 
be a higher voltage transmission network looking down into a more com- 
plex lower voltage distribution network. XS may also be a relat.lvely 
extensive region where power and voltage measurements are not made which 
is siu'rounded by CS where measurements are made. To aid in analysis 
notation OS is further subdivided into "Internal System” or ]S which 
has no immediate connections to XS and "Boundary System” which has 
imm.ediate connections to XS. V?hat is to be determined is an equivalent 
model for XS named "Equivalent System” or ES which will affect CS the 
same as XS. 

It will be assum.ed that measurements are available in OS but not from 
XS. There are two aspects which set this particular identification prob- 
lem apart from many others. One is that the identification must be 
primarily passive. It will not be possible to make any great manipula- 
tion of input signals (power flows) that unduly upset the system or do 
not conform to power consum.ption requirements. Also it can be expected 
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Figure 1.1 Separation of the Entire Netv/ork 
Into Regions 



- 13 - 



that power requirements in XS will be about the same as in OS, and 
since the unknown inputs in XS will cause unknown changes in the 
outputs of CS , "signal to noise ratios" will be about one to one. 
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l. 3 Contents of Thesis 

The basic conventional model which is used for analyzing networks 
will be developed in Section 2 . Variables of interest for this model 
are real power, reactive power, voltage magnitude, and the relative 
voltage phase angle. The m.odel is nonlinear and involves cross coupling 
of all variables. However it will be shown that under certain conditions 
the effect of voltage m.agnitude changes on real power is much less than 
the effect due to phase angle changes, so that an approximate linear 
model can be developed which relates real power to voltage phase angle. 
The linear model is not conventional but it is not new either. One of 
its chief advantages is a very significant reduction of computation 
requirements. But even though a linear model will be used for identi- 
fying the model for Equivalent System, it is still possible to incor- 
porate the identified equivalent linear model in the nonlinear model fcr 
Own System. In Section 3 the linear model for External System is 
reduced to an equivalent model. Then two m.odels for use in identifying 
Equivalent System are developed using the equivalent model. One of the 
models is attractive from the point of view that it only requires obser- 
vations at the tie lines joining CS to XS, but it will turn out that 
there are inherent problems with this model due to the correlation 
am;ong variables at the tie lines. Both m.odels will be placed in a 
general form and in Section 4 the solution to the general problem will 
be found subject to certain conditions. The solution, which has a very 
simple form, is derived from, both an assimied mathematical probability 

m. odel for the disturbances and the method of weighted least squares. 

Error analysis equations are also derived to study the sorirce of errors 
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and find a measure of confidence in the identified parameters. 

Applying the solution of the general problem to the two specific 
identification models in Section 5» the effects of the conditions 
placed on the solution are analyzed. One shortcoming of the more 
easily implemented model is disciiissed in terms of independent inputs. 

A simulation made to tost the proposed methods of identification is 
discussed in Section 6, and the results show excellent correspondence 
with the theory. Errors in the identified parameters are within pre- 
dicted accuracy. Finally, in Section 7 several aspects of the specific 
problem are discussed which involve increasing identification accuracy 
and verification of the identified m.odel. It is also pointed out how 
models for identifying the equivalent system can be valuable by them- 
selves for the purpose of network model reduction. 
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2 . E EECTR IC A L PCTaER TP^ fS M ISS 10 N MODE IS 
2 .1 Ibnlinear Load Flow Rodel 

In reality an electrical power system consists of a large variety 
of devices, many capable of storing energy with time constants and 
natiiral frequencies spanning a very wide range. However for the pur- 
pose of computing power flows in sinusoidal steady state only a few 
types of com.ponents need be considered. The most important of these are 
the transmission lines which transmit power from points of generation 
to loads, the transformers, and the buses which form connecting points 
or nodes for transmission lines, tr a r.s formers , generating stations, and 
loads. The values of interest are the power flows and bus voltages. 
Although the transformers, transmission lines, loads, and sources are 
capable of storing energy, when the system is in sinusoidal steady state, 
the usual approach is to treat energy storage devices as complex 
impedances or admittances and the sources and loads as injections of 
complex power. Hence although the system is oscillating, it can be 
analyzed as a static network. 

A transmission line is of course a distributed param.eter device. 
However as with most such devices, it is possible to model it as a set 
of lumped elements provided the physical size is small compared to the 
wavelength of the power transmitted. Since the wavelength of 60 hertz 
alternating current is roughly the distance between few York and San 
Francisco, while most transmission lines are far less, lumped element 
models are quite valid for this purpose. The pi model, common for this 
type of distributed parameter device, is used here and consists of a 
resistance and inductanice in series with shunt capacitance at each end. 
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Figure 2,1 Pi Transmission Line Model 




Figure 2.2 Electrical Po\ver Network Lumped Model 
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Since the system st\K3ied is assumed to be in sinusoidal steady 
state, the capacitances and inductance can be represented as admit- 
tances. A bar above a symbol vill denote a complex quantity (e.g. 

J^i 

e = E^ cos + jE^^ sin D^). An underlined symbol will 
denote a matrix or vector. Let 



ik 



= y 



-j^ik 



ik 



^ik 



+ jwL 



'ik 



be the complex admittance of the series resistance and inductance, and 
let 



ys 



ik 




be the admittance of the shunt capacitance. The current entering the 
transmission line at the i^^ node is IL^ = E^ ^ ^^i ” ^k^ ^ik * 

Complex power is P + jQ = E I (I is the com.plex conjugate of I) so 
the power enterirg the transmission line from the i^ node is 



’’hk - = ^i"^ik ’ y=ik - ^k’ ^ik 






( 2 . 1 ) 



The network mo(3el consists of admittances which represent the 
transmission lines and transformers connected to buses or nodes. Also 
connected to each node is a current source which represents the load 
or generator at the bus. Power Injections into the buses are positive 
for generators and negative for loads. As shown in Appendix A, a bus 
admittance matrix can be formed which represents the relation between 
the current sources at each bus to the bus voltages. Let Q 

vector whose elements are the complex bus voltages, ® vector of 

the corresponding values of the complex bus current sources, and 
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the coirplex bus admittance matrix. Then the relation is 

--- th - N__ 

L = Y, E, or for the i bus, I. = ^ Y..E, where there are 

-=^us -bus -bus i ik k 

trll 

N buses. Correspondingly, the com*plex power injected into the i bus 
by the current source at the bus is 



P. 

1 



- jQ. = 



E .*2 

^ k=l 



N 



^ ^ •'“k „ -»ik 

2- E^e Y^j^e 



k=l 



- Teey 



k=l 



( 2 . 2 ) 



The load flow problem is that of given bus power injections, find 
the bus voltages and the transmission line power flows. It should be 
noted that even though the network is modeled with linear elements, the 
power flow equations are nonlinear. 
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2 ,? Linear Real Power - Voltage Angle Model 

Under certain assumptions it is possible to reduce the complexity 
of the model. From Equation (2.1) line real pover flow is 



= Real 






-j0 



ik 



i •'ik" 



- E.E 



iVik^ 






= E. y.j^ cos cos(0.j^ + D. - D^) 

^^ik = ^i^ik^°^ ^kpi - \ ®iVik""” ^ik - V 

(2.3) 

If it can be assumed that voltage magnitudes do not vary much from 

nominal operating points, that the transmission line reactance is much 

greater than the resistance so that ^ 90®, and that voltage angle 

differences (D. - D, ) are small, then when E, = (E,) . ., and 

1 k ’ 1 1 nominal’ 

E = (E ) . , PL., E E y.,sin(D. - B, ). For convenience it 

2 2 nominal’ ik 1 2 ik ' i k 

*%» ^ 

will be assumed that E^ = E^ = 1.0. Using the approxim.ation for the 
sine, sin(D^ - D^^) ^ small 

(2.4) 



PL., 

ik 



yik'"! - ®k’ 



As shown in Appendix for such a case a linear load flow equa- 
tion, = -bus“bus’ ^®Lates real power injections to voltage angles. 

For a system with (K+l) busea,. is a Kxl vectoi' of bus real power 

injections, B,^^^ is a Kxl vector of bus voltage angles, and is a 

KxK matrix whose elements are nearly the samie as the magnitude of 
elements of the complex matrix (K+l) bus is the reference 

bus at which the voltage angle is specified. Power injection at the 
reference bus is such that the algebraic sum of real power injections 
of all buses is zero. Rom [ 6 ^ and Rom an^A Schweppe C 8 3 proposed the 
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use of this linear load flow model for on-line estimation of voltage 
angles and transmission line power flows while Baughman and Schweppe C2 ] 
showed how its results compare with the nonlinear model. 



2,3 Discussion of the Linear Load Flow Model 

Although the linear load flow m.odel does not give the exact values 
of voltage angles obtainable by the nonlinear model nor does it account 
for voltage magnitude variations, there are some great advantages to 
Using it for real power flow analyses even when the conditions on low 
line resistance and small voltage m.agnitude spread are only mildly met. 
Foremost is the simplicity of calculation. For (K+l) buses the non- 
linear model requires an iterative solution of a set of 2K simultaneous 
nonlinear equations. The linear model involves K simultaneous linear 
equations. If the tewton-Raphson method is used to solve the first, 
som.e comiparison between the two can be made. The ffewton-Raphson 
method for one iteration is of the fom Aa^ = , Since calcula- 

tions for solutions of simultaiieous linear equations are about propor- 
tional to the square of the number of equations, one iteration for the 
nonlinear miodel requires about four times the calculations for the 
linear model. Also the b'^ matrix must be calculated at each iteration 
whereas linear model remains fixed for a given system. 

Finally, depending on the accuracy required, the Jbwton-Raphson method 
requires about four or m.ore iteratior^. Another advantage is the ease 
with which the linear equations can bo manipulated. This will become 
very important in transforming the network for finding an equivalent 
system as it results in much insight. Historically, before digital 
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corr.puters were in wide use, load flow analyses were ir.ade with calculating 
boards which used physical con'ponents to model the network. The more 
versatile ^.-C calculating boards could account for voltage magnitudes 
and line power losses but were expensive and there were only about 50 
in the United States till The linear load flow m.odel is somewhat 
analogous to the D-C calculating boards which w’ere m.ore numerous and 
less expensive. The linear load flow model accounts for individual bus 
real power injections while some D-C calculating boards used only one 
source . 

The error involved in using the linear load flow model of course 
depends on how well the three conditions are m.et. The assumption of 
small voltage an^le differences between the ends of a transmission line 
usually results in the least error. Except in rather rare circtcr.starices , 
angle differences are within 30° a*^d usually they ere within 10°. Com- 
paring the values of the sine of an angle and the angle: 



1 (degrees) 


Angle (radian-s) 


Sine 


10° 


fg = 0.1745 


0.1736 


20° 


f- = 0.349 


0.342 


30° 


f- = 0.524 


0.500 



Even up to 30° the approximation is within 5^» and under 10° it is 
within 0.5^. 

The error involved in the other two assumptions can be analyzed 
through the nonlinear equation for line power flow. If the line 
conductance, g^^^ = y^^^ cos the line suseptance, = “^ik *^ik’ 

are used then equation (2.3) is 






c 






I 




(2.5) 
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* ^k ^i^k^ik " \) 

In or<^er to make a rather rough estimate of errors let "a” represent 
either or y^i^. Then when is approximated by b^^(D^ - D^) or 



ik* 



ik 



ik' i 



yik(r>i - \), for small angle differences 



■hk ” * ^i'^k*>lk<°i - “k> 



and the relative error for PL., a(D. - D, ) is 

ik ' 1 k' 

11- PL - PL ^ 

relative error = — = c 



- V - E,g,„(E, - E^) - E,E^b,^(D, - 



®l'^k’=ik<'’l - °k> 



a - E ,E, b., 
1 k ik 

1 k ik 



1 g., E. - E, 

^ik 1 k 

\ '’ik 



For E. R* 1,0 and E, r* 1.0, €. is approximately 

X iC 

f ^ y -p ®ik ^i " ^k 

^ ^ vl — - E.E, - r — f; ?7" 

’’ik " ”ik °i - ®k 

If the voltage drop is expressed as a percentage, AE = 1C0(E^ - Ej^)> 
and if angle difference is expressed in degrees, AD = 57*3 (D. " ^i.) 
^ 50(D^ - > the relative error is approximately 



- Z' ® ^ 1 ^ik 



(2.6) 



The purpose of equation (2.6) is to show how the greatest portion of the 
error enters the linear transniission line power flow equation (2.4) • 

g 

The term r — - E E, shows how voltage magnitude changes from nom.inal 
ik ^ 
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values affect the error, while 



1 

2 




&E 

AD 



shows the effect of resis- 



tance anti voltage drop. Table 2.1 compares line flows as computed 
both by the nonlinear and linear equations along with actual errors 
and approximate errors predicted by equation (2.6). It is not immedi- 
ately obvious whether b^j^ or y^^j^ (actually where 

and Ej^ are nominal values) should be used for a in ^ a(D^ - Dj^)* 

In view of the relatively small error and much less involved 
calculatior.s required, and to preserve insight, the linear load flow 
model is used for the purpose of identifying the equivalent system model. 
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3. MODE IS FOR IDE IT IF IC AT ION OF THE EQUr/AIENT SYSTEM 



3.1 Tie Line Pov’er Flow Model 

As previously mentioned, tvo approaches to finding the equivalent 
system model will be used. The tie line power flow model views the 
external system as a separate system and attempts to find an equivalent 
model by usir.g the tie line power flows and biis voltage angles where the 
boundary around CS crosses the tie lines connecting CS to XS . The bus 
injections and corresponding bus voltage angles can be separated into 
3 groups, those belonging to IS (internal System.), those belonging to 
BS (Boundary System), and those belonging to XS (External System). 
Collectively IS and BS form CS (Cwn System). ^ and D_I are the bus 
real power injection vector and the bus voltage angle vector respectively 
whose elements belong to IS. Similarly PB and DB belong to BS , and M 
and DX belong to XS. The letters I, B, and X will usually refer 
values belonging to IS, BS , and XS respectively. Since by earlier 
definition IS and XS are not immediately connected to each other by 
trar.snission lines, the linear load flow equations can be witten: 



”pl~ 




~YII YB 0 




DI 


P3 


= 


YB» YBB Y3X 




ra 






0 YBX ' YXX 




DX 










_ 



(The transpose of a m:atrix will be denoted by Y' = transpose of Y). 
The zero submatrices in the corners are due to the lack of direct 
coupling between K arxl XS. The only coupling between CS and XS is 
reflected by the submatrix YBB . Since the elem.ents of YB3 are values 
of admittances of transmission lines vjhich have at least one of their 
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ends connected to buses in BS, YBB can be separated into that which 
belongs to OS, YBOS . and that which belongs to XS, YTL . so that 
ITB = YBOS + YTL . YTL vrill be a diagonal matrix whose elements will be 
the admittances of the tie lines connecting OS to XS. Equation (3,l) 
can be written 



pf 




'yii yib 0~ 




'of 




lo 

lo 

lo 






FB 
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YIP* YBOS 0 




DB 




0 YTL ^ 




DB 


PX 




0 0 0 




DX 




0 YBX’ Y^ 




DX 


L J 












_ 




^ _ 



'pi' 





'yii 


YIB ' 
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" 0 


P3 




YIB’ 


YBOS 




PB 




YTL DB + YBX DX 


- - 




- 






- 




— _ 



(3.2) 



EC = m'D3 + Y^ DX ( 3 . 3 ) 

Let PTL = XIL D3 + Y^ DX ( 3 . 4 ) 

Tie line pov;er flow, PTL, is a vector of length equal to ra and repre- 
sents the power vjhich flows out of the buses belonging to BS into the 
transmission lines which connect 3S to XS. Solving for DX in equation 
(3.3) substituting into equation (3.4) 



DX = -Y>:X~^YBX'DB + YXX'^PX 

PTL = YTL ra + Y3X(-YXX~^YBX'D3 + YXX~^PX) 

= ( YTL - YBX Y:CX“^X « )P3 + YBX YXX~^PX 



Define an equivalent external system bus adm.ittance matrix, YEOX, and an 
equivalent external bus power injection vector, FEQX . such that 
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YEO;( = - YBX m~^x« 

FEOX = AQX PX 
AQX = - YBX YXX~^ 

Then equations (3.2) and (3.4) can be rewritten 



'll 




~YII 


YB ~ 




~UI 


PB - PTL 




YIB* 


YBOS 




D3 


- « 










— 



PTL = (YTL + YE0X )D3 - FEQX 



(3.5) 

(3.6) 

(3.7) 



(3.8) 



(3.9) 



Equation (3.8) is simply the linear load flow model for CS. By including 
FTL with PB the problem can be solved as usual. Equation (3.9) for FTL 
shows two contributions. The first, ( YTL + YEQX ) DB is the part of FTL 
due solely to the values of voltage angles existing at the buses of BS. 
The second, FEQX is the equivalent power injections in XS as seen by 
BS. It is the bus power injections from ^ that are routed to BS by the 
structure of XS. It would be desirable to know both YEQX and FEQX . but 
even knowledge of just YEQX would be valuable. For example, if the 
voltage angles and line power flows in OS were wanted after a change in 
bus power injections in OS, then if PX (and thus PEQX ) did not change 
between times t^^ and t^, the change could be found by substituting 
equation (3.9) into (3.8) so that at any time t 



'Pj(t) 




"yii 


YB 




"Dl(t) 


PB(t) + FFQX(t) 




YB* 


(BB + YEQX) 




DB(t) 



and for the change with PS'QX (t^) = PEQX(t^) 
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’pl(tp)- PI(t^)‘ 




yn IB 




'51 (tp 


- Dl(t^)~ 


P3(tp)- ra(tp) 




YIB» (YBB + YEQX) 




D3(tj) 


- m(t3^) 



vhere P_T(t^) and 5^(t^) are found by measuring PTL(t^) and using 
equation (3.8). On the other hand, if the voltage angles and line 
power flows in OS were wanted after a change in the status of one or 
more transmission lines in OS, then if YII . YIB . and YB3 are the admit- 
tance matrices after the change, for Pl(t 2 ) = Pl(t^), raCt^) = PB (t . 
and FSQX (t^^ = FEQX (t^^) 



~YII 


YB 




’Dl(t^)“ 




YII YB 




'ei(tj) 


YIB» 


YBE+YEOX 




E3(t;^) 




YB» YM+YEQX 




5S(tj) 



Everything is known except ^^(t^) and ^(tp) since again Dl(t^) and 
DB(t^) can be found by measuring PTL (tj^) and using equation (3.8). 

As formulated, in equation (3.9) BB is an input vector and FTL is 
an output vector while PEQX m.ight be thought of as an unknown distur- 
bance. If a series of input, output m.easurements were taken, then it 
might be possible to find YEQX . Equation (3.9) will be referred to as 



the tie line power flow model. 








* 
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3.2 Boundary Bus linpet^ance i-'odel 

The second approach involves an impedance matrix instead of an 
admittance matrix. In order to obtain the required form, an equivalent 
internal system, as seen by BS looking into B vill be found which is 
similar to the equivalent external system. The top row of equation (3.8) 
is 

= YII ^ + YJB BB 

Solving for ^ 

DI = -yir^YIB ra + YII'^PI 
and substituting into the secoix’. row of equation ( 3 * 8 ) 

P3 - PTI = (-Yir^YB B3 + YII~^PI) + YBCS DB 
Using equation (3.9) for FTL and combinir.g terms 

PB + PEQX - YIB'Yir^PI = ( YBCS + YTL + YEQX - Y 13 ' Y IE )DB (3.10) 



Define an equivalent internal admittance miatrix, YEQI . and an equi- 
valent internal bus power injection vector, PEO I. sim.ilar to ~!:'E9.X and 
PSQX by 



YEQI = -YIB'YII"^YB 


(3.11) 


PEQI = AQI PI 


( 3 . 12 ) 


AQI = -YlB’Yir^ 


( 3 . 13 ) 


Then equation ( 3 .IO) can be rewitten. 

PEQI + PB + PEQX = (YEQI + YBB + YEQX)DB 


(3.14) 



I 
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PEQI is the power injection of ^ which is routed to the buses of BS 
by the structure of B. YEQI is the structure of IS as seen by BS. 

YB3 was previously separated into that belonging to CS , YBC6 and that 
belonging to XS , YTL . If YECS is further separated into that belonging 
to IS, YBIS arj3, that belonging to BS , YBBS so that YB06 = YBIS + YBBS 
or YBB = Y3 IS + YBBS + YTL . then the quantity 

(YBB + YE2I)D3 - FEQI 

is the pow’er flowing out of the buses of BS into the transmission lines 
connecting IS and 3S. If BS had no transmission lines unique to itself, 
then YBBS would be zero. By rewriting equation ( 3 . 14) with 
YBB = YFIS + YBBS + YTL 

PEQI + re + PSQX •= ( YBIS + YEQI )DB + Y BBS ^ + (HL + YEQX )PB (3.14a) 

Equation (3.14a) can be considered to be two matrix Thevenin or Norton 
equivalent circuits, one for IS and one for XS, coupled to the circuit 
for BS as illustrated in Figure 3.1. If equation (3.14) is solved for 
BB, then 

B3 = ( YEQI + YBB + YEQX )~^( PEQI + re + PE OX) 

As shown in Appendix B, ( YEQI + YBB + YEQX ) ^ is that part of the Whole 
System (WS) impedance matrix which relates re to DB. Therefore denoting 
ZBB for that part of the impedarjce matrix 

ra = ZBB ( PEQI + re) + ^ PEQX ( 3 . 15 ) 

In this formulation (P SQI + re) is an input vector, BB is an output 
vector, and ZBB PEQX is an unknown disturbance vector. Again if a 
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RLL=i> 



DB 




FBQI f ^ + PEQX = ( YBQI YBIS )DB f Y3BS D3--t-(YEQX + Y?L)DB 
Figure 3«1 Equivalent EetY/ork Viewed From BS 
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Figure 3.2 Tie Line Power Flow Fodel 
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Boundary Bus Impedance Yodel 
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series of input,, output measurements were taken, it m.ight be possible 
to find ZBB » 

For comparison purposes and in order to be able to investigate 
sources of approximation error later, the analogous equations for the 
nonlinear model are developed in Appendix C. 



3»3 Separetinm Ojn System Estimation From Equivalent System Identification 
The tie line power flow model and the boundary bus impedarice model 
are illustrated in Figures 3.2 and 3.3. As shown in Table 3.1> both 
approaches result in equations of the form z = H u + y where u is the 
input vector, z the output vector, y a vector of unknown disturbances, 
and H the structure m.atrix to be identified. Both z and u were assum^ed 
to be known. However, since the state of OS (i.e. the bus voltages) can 
be estiiTjated by knowing what power enters from XS thi'ough the tie lines, 
the problem of estim,ating the state of OS can be separated from, identi- 
fying the parameters of the equivalent system. What happer^s to CS 
depends on XS , but the effect of XS is sxmm,arized by the power flowing 
through the tie lines. Hence existing state estimation techniques C 5 ), 

£8 and [9 1 can be applied to obtain z and u. In practice there will 
be soir.e error associated with metering pov/er and estimating voltage 
angles. In such a case let v be the error of outout estimate and v 
be the error of input estimate. Then 

(z + V ) = K(u + V ) -f V 
- -2 - - -ni - 

or z=Hu+(v-v-fHv) 

- ~2 U 
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Table 3.1 

Identification Kodels 
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By replacing v - + H with a new unknown disturbance vector y 

z = H u + y 

the same approaches can be used. The only difference is that there is 
generally more disturbance associated with identifying H. 



3 . A Use of Changes in Variables 

Up to now z, u, and y have not been associated with time. From 
now on z(t), u(t), and y(t) will be used to designate a set of z, u, 
and y for a particular time t. The disturbance vector y(t) will be 
treated as a noise term since u(t) is analogous to an input signal and 
z(t) is analogous to an output signal corrupted by noise. Ordinarily 
when working with estimation or identification problems, it is much 
easier if noise termc have a zero mean. In the two approaches above 
neither will in general have zero temporal m.ean noise vectors v(t). It 
is expected that y(t) will vary about some value depending on the time 
of day, the season, etc. However, if the differences in z(t), u(t), 
and y(t) are taken between two different times, then it may turn out 
that the difference or charge y(n) = " -^^n^ will have a zero 

mean. The general equation is then 



y(n) = H u(n) + y(n) (3-16) 

z(n) = 
u(n) = 
y(n) = 



where 
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It will turn out that there are other reasons for using the difference 
in variables rather than the actual values of the variables. The most 
significant is that using the actual values may result in subtracting 
variables whose magnitudes are very large compared to their difference. 

From now on, any reference to the previously defined sjnnbols for 
voltage angles or power will refer to the changes of variable between 
two times rather than the actual values of the variables. Hence the tie 
line power flow model is 

FTL(n) = (YTL + reOX)DB(n) - FEQX (n) (3.17) 

and the boundary bus impedance model is 

^(n) = ZB3 [ PEOl (n) + ra(n)] + ^ FEQX (n) (3.18) 

where PTL(n) = FTL(t^^) - FTL(t^) 

ra(n) = m(tj^j_) - ra(t^) 

53(n) = ra(t^^) - DD(t^) 

FEQl (n) = KOl(tj^^) - PEQl (t^^) 

PEQX(n) = PEQX(t ,) - FEQX(t ) 

^ ^ nfl n 

for N sets of measurements n = 1,2,...,N. 
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4. IDEIiriFICATION CF LINEAR, STATIC SYSTE>S 
4.1 Conditions on the System 

Both formulations for identifying the structure of the equivalent 
system have resulted in the same general equation, z(n) = H u(n) + v(n) 
where there are N sets of input vectors, u(n), and output vectors, z(n). 

H is the structure matrix and v(n) is a pseudo noise vector. The 
general problem, can be stated: 

Given: A set of input ar.d output vectors u(n) and z(n), n = 1,2,...,N 

Find: The structure matrix H and statistical information on y(n). 

Several assumptions on u(n), v(n), and H will be made and later 
these will be justified when the results are applied to the specific 
problem;. First elem.ents of H are assum.ed to be independent. Specifically 
it is assvim.ed that H is not symmetric. This is certainly not true for 
the specific problem, but various justifications will be given later. 
Second, y(n) is assum;ed to be a zero temiporal mean process, or 

= 0. where E^a| is the expected value of a. Third, v(n) is 
uncorrelated in time, or E^y(n) y'(m.)^ = 0 for n / m (n represents a 
sequence in time). Fourth, the covariance m.atrix for v(n) is assumed 
to be of the form R, or E^(n)y' (n^ = R(n) = 

be called a confidence coefficient for the n measurement set and 
indicates the relative noise level for that set. Larger c(n)^s indicate 
less noise. Fifth, v(n) is uncorrelated with u(n), or E|x(n)u’ (n)^ = 0. 

At this point there are many possible approaches to the problem. 

Two will be used here. One is a time tested mathematical approach 
which assumes a particular probability model for the noise, and the 
other is an engineering approach, also time tested, which uses no model 
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for the probability distribution of the noise. It will turn out that 
both result in the same solution which will have a relatively simple 
form. The two approaches ere: 

1. Assume that the noise, v(n), has a Gaussian probability 
distribution (with a zero mean and covariance of - 7 — c R). 

Then find the miaximum likelihood estimates of H and R. 

2. Assume no particular probability model for the noise and 

A 

use the miethod of weighted least squares to find H. Then 
m.ake a reasonable estimate for R. 

The fact that the two approaches result in identical solutions is a 
well known corisequence of using quadratic criteria for optimization 
and estimation of linear systems. However, in addition to analj't.ical 
convenience, it is also well known that in m.any instances the Gaussian 
probability distribution itself is a rather reasonable probability 
model. 



4 .2 yaximvm Likelihood Identification 

First the solution will be found for the m.axim.um. likelihood approach 

assuming Gaussian noise. There are four known items: the set of input 

vectors u(n), n = 1,2,. ..,N, the set of output vectors, z(n), n = 1,2,. ..,N, 

the correspoTiding confiderce coefficients, c(n), n = 1 , 2 ,...,N, and the 

form of the probability distribution function for v(n), n = 1,2,...,N, 

which is zero mean Gaussian, with a covariance matrix 7 R. Note 

c^nj - 

that R itself is not known. Therefore if estimates of H and R can be 
found, the probiem is solved since knowledge of R com.pletely specifies 
the statistics of v(n) when y(n) is Gaussian. 
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The probability distribution function for the Kocl vector v(n) is 

r 1 r K -K 1 “p -p X’ (n)c(n)R“^v(n) 

p[x(n)] = [(?"■) c (n)lR)J e ^ 



The likelihood function of z(n) for a given H and R when z(n) and u(n) 
are Kxl vectors and H and R are KxK matrices is 

_1 

p[z(n),u(n):H,R] = [(2rr) c“^(n) \Rl') 

where 

j[z(n),u(n) :H,R] = |[z(n) - H u(nj) c(n) R~^[z(n) - H u(n)] 

Since v(n) is uncorrelated in time (E^v(n)v' (mj^ = 0 for n / m) , the 
joint probability distribution function for the set V = [v(l),v(2),...,v( 
is the product of the distribution functions for v(l) ,v(2). . . ,y(K). 



p(l) = f[x(1)] p[x( 2)] ... p[v(N)] 



1 1 

n=l 



2 v’ (n)c (n)R~\(n) 
n=l 



Let p(Z:K,R) be the likelihood function of the sets 
Z = ^z(l),z(2),... ,z(N)]-, and U = ^u(l) ,u(2) , . . . ,u(N)| for a given H 
and R. Then 

_1 

p(Z,l!:H,R) = [(2.)® IR1« H =-"(n)]’' 3-^(S.a«,S) 

n=l 

where 

J(Z,U:H,R) = ^ ZJ[z(n) - H u(n)] c(n) R‘ [z(n) - H u(n)] 
n=l 

A A 

Let H and R be the values of H and R v/hich max3.mize the likelihood 
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function when Z and U are the actual measured sets. As developed in 
Appendix D these maxirnum likelihood estimates are 



A 

H = 



2] c (n)z(n)u' (n) 
n=l 



X 1-1 

X c(n)u(n)u' (r.)j 
n=l J 



(4.2) 



and 



N 

R = ^ X c(n)rz(n) - H u(n)][z(n) - H u(n)] 
n=l 



(4.3) 



provided the vectors u(n), n = 1,2,...,N are such that the inverse of 

r N 1 

12 c (n)u(n)u* (n)l exists. 

ln=l 

Two things should be noted. First, it is not necessary to know 



A A 

R in order to find H. This is significant in that it is not necessary 

A A 

to solve for R and H simultaneously which would most likely result in 

A 

greatly increased computation requirem.e nts . The second is how R varies 
with H. The quantity K|trace (R)] is actually 2J(Z,U;H,R) when R is set 

A 

equal to the identity m.atrix. As seen in Appendix D, H is the value of 
H such that the m.atrix gradient of J with respect to H is zero. Then 



Itrace(R)^ = ]2J ] = 2 



or in some respect R is not sensitive to sm.all errors in H. 












I 

It 
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4>3 Least Squares Identification 

The Gaussian distribution of v(n) was assuir.ed in order to provide 
a probalistic framework, however the sam.e equations can be obtained by 
assuming no probability model for v(n) and instead relying on the 
nearly universal method of weighted least squares* The sum of weighted 
square errors between z(n) and K u(n) is 



where is an unknown weighting matrix. Since S is the same as 



estimate for H did not depend on R any\^ay and the value of H which 
minimizes S for given sets z (l) ,z (2) , . . . ,z ( N) and u(l) ,u(2) , . . . ,u(N) 
will be the sajne as the maximum likelihood estimate for H. Under such 



S = I 2 - H u(n)] c(n) T [z(n) - H u(n)] 



n=l 



J(Z,U;H,R) when T = R from Appendix B it can be seen that the 



circimstances , a reasonable 






n=l 



which is the same as the maxiir.vun likelihood estimate equation (4.3)* 







I 
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4.4 Error Analysis 

The error analysis which follows only uses the covariance of y(n). 

For both the approach assming Gaussian noise an-d the least squares 
approach, the covariance of v(n) was assrimed to be E^v(n)y • (n)^ = R. 

Therefore results of the error analysis are applicable to both approaches. 
W.en the equation for z(n) is substituted into equation (4.2) 

A r N N -vr N 1-1 

H = I 2. c(n)H u(n)u'(n) + ^ c (n)y (n)u* (n) 11 ^ c (n)u(n)u» (n)l 
Ln=l n=l ]ln=l J 



N 

M “ M = 1 2 

ln=l 



N 1 

2 c (n)u(n)u* (n)l 
n=l J 



-1 



(4.4) 



One of the conditions on the solution was that v(n) was not correlated 
with u(n). It is interesting to determine the effect when that condi- 
tion is not met. Let W be a ’.weighted average of u(n)u*(n), n = 1,2,...,N 



W = 



c 

o 




-1 • 


Ln=l 





N 1 

2. c(n)u(n)u» (n) 
n=l J 



(4.5) 



If y(n) is correlated v;ith u(n) so that the weighted average of 
y(n)u*(n), n = 1,2,. ..,N, is 





N 


-1 


N 


L- = 






2 ® (^)z(n)H* (n) 




n=l J 




.n^ J 



^|z(i'‘)y' (n)"^ 



Then 



A _1 

H - H = M W 



(4.6) 

(4.7) 



and as N increases H - H will not go to zero as would be the case if 
y(n) were not correlated with u(n). Error is introduced in the form 
of a bias which is proportional to the correlation of v(n) with u(n) 
and inversely proportional to the "power" of u(n). 

From equation (4.2) it m^ay be noted that only the i element of 
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z(n) affects the estimate of the element of H. Similarly from 

til 

equation (4.4) only the i element of y(n) affects the i row of 

A 

the error H - H. This reflects the fact that K was essuired to be, 

th 

in general, not symmetric, and that the i element of z{n) is affected 

•f V» + Vi ^ 

only by the i row of H and the i element of v(n). Let h^ be the 
i^^ row of H. 

H = 




Then the estimate of the 



.th 

1 



row of H is 



A I 



h. 

-1 



■ N 

^ c(n)z (n)u’ (n) 
n=l 



' N -| -1 

2 c(n)u(n)u» (n)l 

[n=l J 



(4.8) 



arid the error is 



\ t I 

h. - h. = 

-1 -1 



■ N 


■ N 


^ c(n)v. (n)u» (n) 


V c (n)u(n)u' (n) 


n=l ^ 


,n^ 



-1 



(4.9) 



th *th ^ 

The error covariance between elements of the i row and j row of H 



is 




I 
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= E 



^ ^ c (n)u(n)u* (n)| c (n)u(n)v^(n)| 

^ c(n)vj (n)u' (n) c (n)u(n)u» (n)j 

N -1-1 r N N 1 

c(n)u(n)u» (n) 1 ^ 2 c (i’’)c (n-)ECv. (n)v . (m)|l 

n=l J Ln=l in=l i i J 



■ N -j -1 

X 2 c(n)u(n)u' (n) 

Ln=l J 

Since E^v(n)v’ (m)^ = 0 for n / m and E^v(n)v' (n)^ = R 

- llj) ^ 1^2 c(n)u(n)u» (n)^ 






( 4 . 10 ) 



Let c be the average of the c(n) ’s. 

1 

c = ■», S c(n) 

n=l 

Then in terms of W, the weighted average of u(n)u* (n) defined by equation 
(4.5) the error covariance is 



P. . 
“ij 



1 

N 



R. . , 

C - 



(4.11) 



As might be expected, on the average, the error covariance of the esti- 
mate is inversely proportional to the number of measurement sets. In 
terms of signal to noise ratios, the quantity c ^R. . is a mieasure of 
the power of the noise v(n), and IH is a m.easure of the power of the 
signal u(n) . 





1 
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5. IDENTIFICATION OF PQ-JER SYSTEIS 



The problem of identifying the equivalent admittance matrix and 
estimating statistics of the equivalent bus power injection vector 
for the external system model has been placed into a general form, 
and the solution to the general form has been given subject to certain 
conditions. The two models for Identifying the equivalent system are 



where n = 1 , 2 , . . . , N. 

In both models variables such as ra(n) are the changes which occur 
between tijmes t^ and M(n) = DB(t^^) - D3(t^). The general 

problem, is 



The Tie Line Power Flow Model 



PTL(n) = (YTL + YEQX)DB(n) - PEQX (n) 



The Boundary Bus Impedance Model 



D3(n) = ZBB [ PEQl (n) + ra(nj\ + ^ PEQX(n) 



z(n) = H u(n) + y(n) n = 1,2 



N 




and the solution 



H 



A 




R 



c(n)^z(n) - H u(n)j [^z(n) - H u(n)] 



n: 



was based on the assuifiptions 
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1. Eler.ents of H are independent of each other. 

2. E^v(n)| = 0 

3. E^v(n) V* (roj^ = 0 for n / m 

4. E^v(n) v’(n)^ = R n = 1,2,...,N 

5. E^v(n) u* (n)^ = 0 

ibw the effects of these assumptions on each model will be discussed. 

For both models if symmetry of H were taken into account at the 
beginning it would be necessary to perform matrix inversions of dimen- 
sion ^ K(K+l) where K is the dimension of H. In developing the error 

A 

covariance equation for K it was noted that with H not symmetric , each 

row of H is actually identified independently of the other rows. of H. 

With sjTTjretry this is no longer true and the result is that all irde- 

pendent elements of H (either upper or lower diagonal must be identified 

together. The size of the matrix inversion necessary can be significant 

as the dimension of H increases. For example, if the dimension of H 

is 10, the dim.ension of the matrix inversion is 10 when H is not 

assumed to be s^/mmetric because as was shown earlier, each row of K is 

identified independently of the other rows and only one matrix is 

inverted. However, an inversion of dim.ension 55 is required when 

symmetry is accounted for since the 55 independent elements m.ust be 

identified together. One way to handle symmetry at the end is to 

treat the ij^^ element and the ji^^ elem.ent of H as two estimates of 

the ij^^ element of H (for i / j). Then using the corresponding error 

til 

variances, the weighted estimates can be combined. Let a be the ij 

^ th ^ 

element of H and r its error variance, and b the ji element of H 




m 
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and q its error variance, then the combined estimate c is 



c = 



a 

r 



1 

r 



+ 



+ 



b 

q 

1 

q 



SjTranetry of ZB3 and YEQX is merely the result of reciprocity since 
they represent passive networks (the sources, FEQI . P3 . and PEQX have 
been moved outside). However, in the esse of ZBB what is ultimately 
important is the product ZBB ( FEQI + P3). By the principle of recip- 
rocity an input at terminal 1 causes an output at terminal 2 equal to the 
output at termiinal 1 caused by the same input at terminal 2. But because 
of different dem.ands throughout the network, the input at terminal 2 may 
never be as large os the input at terminal 1. 

Had the actual power injection and voltage angles been used, z('tfj) 

(i.e. PEQX(t ) or ZBB PEQX(t )) would most likely have a non zero mean 
n n 

value. However, because changes in values were used, it may be expected 
that v(n) will have a zero mean value provided all the m.easurements are 
not taken when the entire system is m.oving in the same direction. For 
the case when the entire system is moving in the sam.e direction, it may 
turn out that E|^v(n)y ' (n-k)"^ = 0 for k greater than 1 or 2. That is, 
v(n) is only correlated with only the last one or two mieasuremjent sets. 
This could be accounted for by a modification to the approach, but it 
is not anticipated that time correlation of y(n) will be a m.ajor source 
of error. Also, when the entire system is moving in one general direc- 
tion, depending on the sampling interval, there will still probably be 
some independent movement of certain buses in the opposite direction 
due to base loaded generators, load fluctuations, etc. 




i 
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The fonn of R(n) = E^v(n)v* (n)"^ = R was based on the 

assir^ption that the covariance of the noise is dependent only on the 

time interval between measurement sets z(t ), u(t ) and z(t ^), u(t ^). 

D n j. rH'jL 

It is then reasonable to assime that the covariance of v(n) = 
is of the form R(n) = (t , - t )R. There may be error in the measure- 
m.ents z(n) and u(n) but it was shown how such error could be lumped with 
y(n). If y (n) and v (n) are the errors of z(n) and u(n) respectively, 

Z "”U. 

then z(n) = H u(n) + y(n) 



where y(n) = -y (n) + H v (n) + y(n) 

Z U- 

Since y(n) is independent of E^x(^)x' = 2 + cTTiy - 

where 2 = E x^(n) - [m and should be 

relatively constant since it depends on measurement error. The m.easure- 
m.ent error should be in the vicinity of one to ten percent while y(n) 
can be expected to be in the vicinity of ICO^ of H u(n) so the effect of 
2 on the estimate of R should be very minor. 

The last assum.ption, that y(n) is uncorrelated with u(n) can be a 
source of considerable error when this condition is not net. The error 



between H and H was shown to be 



H - K = 



N 



^ c(n)y(n)u' (n) 
Ln=l 



■ N n-1 

2 c(n)u(n)u» (n) \ 
n=l J 



For the boundary bus impedance model 

DB(n) = ZBB[ PEQl (n) + PB(n)] + ^ PEOX(n) 

y(n)u'(n) = ^ PEQ X (n)[ FE0.l (n) + PB(n)V 

= ZB3 AQX[PX(n)PI' (n)A2I' + K(n)ra»(n)] 



while for the tie line power flow model 








I* 
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PTL(n) = (YTL + YEQX )D3(n^ - F£QX (n^ 



v(n)u'(n) = - PEQX (n)D3' (n) 



= RC(n)[^ Pl(n) + ra(n) + ^ ^(n)]* ZBB 

= -AQX [PX(n)PI» (nUQI* + FX(n)ra’ (n)] ^ 

- ^ PX(n)PX« (n)AQX«ZPB 



■where AQX and API were defined in equations (3»7) and (3.13). Both 
models have the term AQX ^PX (n)PT' (n) AQI + PX(n)^’ (n)\ in common. 
However, the tie line power flow model has the additional term 
AOX PX(n)FX' (n)AOX' = PEQX (n) PEOX ' (n). Ordinarily it would be expected 
that a charge in power injection at one particular bus will be indepen- 
dent of the changes in power injections at all other buses so that 
E{pX(n)PJ' (n)^ = 0 and e[pX (n)^ » (n)"^ = 0. The exceptions might occur 
when the entire system is m.oving in the same direction as discussed 
earlier. However in the case of the tie line power flow m;odel, 
E^FEOX(n)PEQX ' (n)"^ = R is not zero unless the external system bus 
injections do not change. This causes a bias in the estimate of H which 
depends on the magnitude of FEQX (n) in relation to FEQl (n) + P3(n), but 
in general it can be expected that the two will be of comiparable m.agni- 
tude so that the bias introduced is in the vicinity of lOC^. Using W 
and K as defined by equations (4.5) and (4.6) 



[ N 3 


1-1 r 


E c(n)’ 




ln=l 


1 l> 


rx 1 


-1 ■ 


- 2 c(n) 




[n=l J 


1 L 
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Z 
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c (n)y (n)u' (n)j 
c(n) FEQX (n) PEQX ' (n) 
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the bies in the estimate for the tie line power flow model is 
H - H = M ZB3 U~^ 

A A 

SO that the estimation equations for H and R are really coupled for 
this model. 

The boundary bus impedance model and the tie line power flow model 
are two specific approaches to the problem of identifying the equivalent 
external system model. In the former the inputs are Own System bus 
power injections and the outputs are boundary bus voltage angles. In 
the latter the inputs are boundary bus voltage angles while outputs are 
tie line power flows (which deperids on voltage angles). Because it is 
power which is bought and sold, while voltage angles are determined by 
the line admittances and the distribution of bus power injections, it 
can be seen that bus power injections are the legitimate independent 
inputs and not voltage angles. This is the cause of the bias in the 
estimate for the tie line power flow model. Both models need D3(n) but 
the tie line power flow model uses the tie line power flow vector, 

PTL(n) ^ which in general will be more accurate a measurement than the 
estimate for the equivalent injections from 13 , FEQl (n) . It is possible 
to m.easure PTL (n) directly while the actual value of PEQl (n) is dependent 
on voltage magnitudes and transmission line resistance as shown by the 

A A 

nonlinear analysis in Appendix C. However the solutior^ for H and R 
are more readily fouad for the boundary bus impedance model. 

Obtaining the estimate of YEO.X from the estimate of ZBB is straight- 
forward. However, if the estimate of YE OX is to ultimately be used in 
a linear load flow FiOdel for CS , then as shown in Appendix B, the 
estimate for ZBB can be used directly. If Z06 is the bus impedance 
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iT'atrix for OS, 







zos 




ZI 


DB 








P3 + FEQX 


— 
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— ^ 
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6 . S 10 N RES ULTS 

In order to test the method proposed, a coiTiputer siiriulation 

program was developed to generate data, identify the equivalent model, 

and check how well the identified model could predict line flow 

charges in OS. Figure 6.1 gives a broad overview of how the simulation 

was conducted. First, using transmission line data, matrices to be 

used in the simulation were computed. These include 2 /a^S (impedance 

matrix for VJhole System), AOI, ¥H~^. ( YEQI + YEffi ) , (YEOX + YTL.) , AQX . 

and YTL. For convenience, the confidence coefficients, c(l), c(2),..., 

c(N), were all unity. Matrices CHI and SGI IV accumulated a running sum 

of z(n)u' (n) and u(n)u' (n) respectively. In order to simulate the 

variations in bus power injections throughout the netv;ork, P(n) (vector 

of bus power injection changes for VS) was generated by random, mjcnbers 

having an PdS (root mean square) value which was a fraction of the 

nominal operating value of bus po’,.'er injections fraction, 

FCTCS, was used for all buses in C6 and another, FCTXS , v;as used for 

all buses in XS. (i.e. E^FX?(n)I = fpCTXS x FX 1^) This main- 

c 1 O ^ nom^J ' 

tained a relative scale for changes while at the same time it simulated 
independent changes in power injections. P-s shown in Section when 
the linear model is used, the bus voltage angle vector for OS can be 
found using either 



Zl(n) 




Pl(n) 


D3(n) 


= Z:1S 


ra(n) 


FX(n) 




PX(n) 



or equivalently 
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Figure 6.1 Basic Plow Diagram For Siimilation 
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Ordinarily equation (6.?) or its nonlinear analog would be used to 
estimate the state of OS using real measurements of Pl(n), ^(n), and 
PTL(nO (and any other redundant r.easurem.ents) . However, for purposes 
of simulation equation (6.1) would have to be solved first to find 
DX(n) in order to generate PTL(n). Since either method results in 
the sarnie values of D3(n) (in the linear model), equation (6.1) was 
used to simulate the output of a state estimation for DB(n). The 
simulated state estimation for OS results in errorless data. However 
as discussed earlier the only effect of mieasurem.ent noise and state 
estimation errors is to modify the disturbance vector v(n). 

Everything described up to now was for the purpose of generating 
data for the identification routine. Using generated data, Pl(n) , 
PB(n), and C3(n), the input and output measurem.ent vectors were formed 
u(n) = AQI Pl(n) + PB(n) and z(n) = DB(n). These values were also 
stored for later use in finding R. Running sums of z(n)u* (n) and 
u(n)u'(n) were maintained by CHI and SGItY so that at any tim.e the 
current estimate H could be found by H = CHI SGIW Using H and the 
stored measurement sets, an estim.ate R was found. Parallel to all of 
this, the actual disturbance vector, v(n), found by v(n) = Z33 FEQX (n) 
= ZBX ^(n), was used to maintain a running sum of y(n)y'(n) in R. 

/V ^ 

The estimated values H anH R could then be compared to the actual 
values H = Z3B and R. In order to compare how well the identified 
m.odel could predict changes in line power flows, the estimate for ZOS . 
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the Impedance matrix relating OS bus power changes to OS bus voltage 
ar^le changes was found by 



zos = 


YIl'^ 0 




491' 


^ r -1 

ZB3 [AQI l] 


(6.3) 




i 

lO 

lo 
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Again random variations of bus power injections were generated using a 
nominal base value and estimated bus voltage angle changes found by 



'A(nf 




'piCn)" 


D3 (n) 


= 


PB(n) 



were used to estimate line flow changes. The actual bus voltage angle 
changes using the same bi::s power changes were found from B(n) = g//S P(n) 
and were used to find actual line flow changes which could be compared 
to the values found using the identified model. 

Throughout the program, in order to present data in a condensed 
form, for mjany variables (z(n), u(n), y(n), actual line flows, errors 
in predicted line flows, and generated bus power changes) only the RNS 
value, the m.ean, the minimimn, and the m.aximum values were maintained. 

It should also be pointed out that the basic approach is simple. Only 
a small portion of the program was involved in the actual identification 

A A 

of H and R. The rest is support for generating data and evaluating 
results. 

Using equation (4.11), the estimated error covariance matrix, 
of row i of H is for c(n) =1.0 
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where 





K X SGIW 



N X S 1GV:A 



Hence the estimated standard deviation of element ij of H is 

(R.. X SIGyA..)'^. Since R. . should be somewhat proportional to the 

11 jj 11 t' t' 

square of PCTXS while SIGKA.. should be inversely proportional to the 

J J 



square of PCTCS, in the case of the linear model with perfect m.easure- 
ments of z(n) and u(n), identification error is dependent on the ratio 
of PCTXS to PCTCS (percentage variation of XS bus powers to percentage 
variation of OS bus powers). 

One system presented in a paper C 3 was used here to test the 
method. The system shown in Figure (6.2) consists of 18 buses. For 
purposes of simulation Whole System consists of these 18 buses of which 
buses 1,2,3 >4, 18 belong to IS , 5 through 7 belong to BS , and 8 

through 17 belong to XS. Bus 18 is the reference bus. Other parti- 
tionings could be used but this one does result in matrices ( YEQX + YTL ) 
and ( YEQ I + YBCS ) which are comparable in size. This avoids a case in 
which couplir.g of BS with IS is very m,uch stronger than BS with XS. 

For this system (YSQI + YSOS ) and ( YEQX + XIL) are 
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Figure 6.2 Ketv/ork Used For Simulation 
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Both the effects of increasing disturbances and increasing 

rceasureirents are simulated here. In Figure 6,3 and Table 6.1 are the 

results of increasing disturbances (PX(n)). The RhS of bus power 

changes for CS in each case was 10^ of nominal while Rl'S values of 

PX(n) were C^, 1.25^, 2.5^^, 5.0^, lOf., and 2C? of Another way 

to say this is that "signal to noise ratios" were oo :1, 8:1, 4^1, 2:1, 

1:1, and 1:2 respectively. The sar.e sequence of random numbers were 

generated each time. Included in Table 6.1 are actual values of ZBB . 

actual errors in identifying it, and the estimated standard deviations 

using R in equation (4.11). Also included are the values of 
, ^ 

( YSQa + YTL ; estimated from ZBB . and the errors of predicting changes in 
line flov;s using ZBB . In a similar format, the results of increasing 
the niimber of measurements are summarized in Figure 6,4 and Table 6.2. 
Changes in both OS and XS bus power injections had R>S values of lOjJ 
of the nominal operating point while H was identified using 16, 32, 64, 
128, 256 , and 512 measurements. 

The results show that not only is identification error commer^urate 
with the noise and num.ber of measuremients , but also estimated standard 
deviation gives an accurate m.easure of the identification error. 

Because the same random, numbers were generated for each value of PCTXS 
in the first case, the actual en*ors shown in Figure 6.3 are nearly 
exactly linear since as was shown earlier, for equal num.bers of 
mieasurements identification error depends on the ratio of PCTXS to 
PCTOS. For the case of increasing mieasurem.ents in Figure 6.4, ZBB^^ 
does not have a steadily decreasing error as with ZHB^^, but it is still 
in the range of the estimiated standard deviation. 
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Ratio OF XS SMS’ CHflWC-ES TO O S BUS CHANC,ES 



Figure 6,5 Identification Error Foi* Increasing 
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To contrast identification models, in Table 6.3 the results of 
identifying YEQX + YTL by the tie line power flow model along with the 
estimated standard deviations are compared with + YTL as found 

from Z33 which was identified by the boundary bus impedance model 
using the seme data. Besides identifying the 'wrong values, the tie 
line power flow model gives inaccurate estiriates of error. 

In Figure 6.5 is another network presented in a paper by Stagg C9 3. 
OS consists of the higher voltage network and XS is the lower voltage 
network. The ’’tie lines" are really tr a TiS formers . Using the bo’andary 
bus impedance model an equivalent network is identified to replace the 
21 buses of XS, so that OS can be analyzed with 9 buses instead of the 
original 30. The results in Table 6.4 again show close correspondence 
between actual error and estimated standard deviation. 

The sim.ulation carried out was certainly not fully realistic. 

However the realistic aspects include 

1. Values of bus power changes were found first, then a state 
estim.ation for OS vjas simulated and appropriate values were 
passed to the identification portion of the program., 

2. The effects of increasing measurem.ents were studied under 
the assimiption that XS would ’wary as much as OS. 

The simulation was unrealistic in that 

1, The nonlinear model was not used to generate data with the 
result that a linear identification model was used to find 
parameters to a linear system. To generate data from a 
nonlinear m.odel properly would require a complete load flow 
for each m.easurem^ent set. (it would be simple to vary voltage 
m.agnitude and voltage angle and then find bus power, but bus 
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power nmst be the specified variable as it is the independent 
input, not voltage angle.) 

2. Bus power injections were independently generated whereas 
realistically they will probably be correlated in space and 
time. There were no conditions on the solution for correla- 
tion in space. In fact FEQX (n) itself has a covariance matrix 
^ ^ PEQX (n) FEQX ' (n)"^ which is in general not diagonal. However, 
the solution did assume no correlation in time. As discussed 
earlier in Section 5, this effect can be accounted for by 
modifying the approach. 

Based on the results of these and other simulations certain con- 
clusions can be made about the boundary bus Impedance model. First, 
with enough measurements it is possible to identify the equivalent 
miodel under conditions of the simulation (independent m*ovement of the 
bus injections). The number of measurements m.ay seem rather large but 
of course the noise is also rather large. The results show that as 
expected either a decrease in the disturbances by a half or using four 
tim>es as m^any measurements approximately halves the error. Second and 
possibly even more im-portant, the error predicted by using the estimated 
R in equation (4.11) agrees very well with the actual answer. This 
mieans that an accurate estimate can be predicted from only the m.easure- 
ments used to identify H so that the limits to which H can be trusted are 
also known. This is of course not the case with the tie line power flow 
model. 

A very important question v/hich m.ay arise is how much data would 
be required to identify the system. The sim.ulation was carried out 
under idealized circumstances and as was shown earlier, for the 
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simulation the ir^entification error for equal nui^bers of measureir.ents 
depenc^s only on the ratio of XS bus power changes to OS bus power 
changes. In actuality however there are upper and lower ILmits, If 
data were used fron samples taken too close together in tim.e so that 
changes vjere small, the effect of measurement noise and state estima- 
tion error could become very significant. On the other hand, the 
system can vary over only so much of a range, so that if the time 
between measurements were too large, then because of cyclical patterns 
of power requirements, there might not be much new information after 
the first few measurements. To determine the amount of data required 
and the best interval between sam.ples would require a study using 
actual operating data. 
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Table 6.1(b) Identification Error For Increasing Disturbances ( YEQX + YTI. ) 
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Table 6.1(c) Mentification Error for Increasing Distiirbances (Line Flows) 
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Table 6.2(a) Identification Error For Increasing Measurements ( ZBB ) 
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Table 6.2(c) Identification Error For Increasing Measurements (Line Flows) 
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YEOX + YTL 



Actual Value 


11.50 


-6.14 


-5.36 


Error BBIM 


.00 


-1.37 


1.30 


Error TLPF]>' 


.30 


-14.44 


10.66 


E.S.D. TLPFM 


.19 


.40 


.39 


Actual Value 


-6.14 


14.26 


-8.12 


Error BBIM 


.07 


-1.87 


2.97 


Error TLPFM 


4.52 


-19.37 


14.18 


E.S.D. TLPFM 


.24 


.52 


.51 


Actual Value 


-5.36 


-8.12 


13.47 


Error BBIM 


.52 


- .52 


2.79 


Error TLPF1-! 


5.41 


-23.24 


16.07 


E.S.D. TLPFM 


.60 


1.29 


1.26 



Actual Value — Actual Value of YEQX + YTL 

Error 3BK -- Error in Identifying ITOX + YTL using the Boundary 

Bus Ihpedance Model YEQX + YTL = ^ - ( YEQI + YBCS) 



Error TLPFH — Error in Identifying YEQX + YTL using the Tie Lin^ 
Power Flow Model 



E.S.l). TLPFM — Estimated Standard Deviation using R from the Tie Line 
Power Flow Model 



OS and XS varied equal amounts (IC^ of nominal) and the same 
256 measurements used for both models. 



Table 6.3 Comparison of Identification Models 



I 
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ZBB for 30 Bus S 2 /stem 



Actual Value 


.C644 


• C459 


.C491 


Actual Error 


-.C€16 


.CCC4 


.G179 


E.S .D. 


.CCI 9 


.CC20 


.0217 


Actual Value 


.C4E9 


• C 665 


.C 663 


Actual Error 


-.COI 9 


.C0C7 


.0184 


E .S .D . 


.CC21 


.GC23 


.0244 


Actual Value 


.C491 


.C 663 


.1122 


Actual Error 


-.C021 


.GOlO 


.0222 


E .S .D . 


.CC25 


.GG26 


.0282 



E.S.D. -- Estimated Standard Deviation using R. 

/\ 

ZBB Identified by the Boundary Bus Impedance Model 



IC^ variation of bus pcvers in both XS and OS. 



256 rieasurer.ent sets used. 



Table 6.4 Identification of a Lover Voltage Ibtwork 
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7. USE Ii: P(7.'/FR SYSTEUS— SOI'E: PRACTICAL ASPECTS 
7 .1 Use of Bus Power Inputs 

From equation (4.1C) which relates the error covariance of the 
estimate, it is obvious that the input directly affects the error. 

For the boundary bus impedance model u(n) = PSQ I ( n^ + PB(n). As 
discussed previously and shown in Appendix C, FE3l (n^ really depends 
on the operating point in a nonlinear manner so that it is possible for 
FEQI (n^ to introduce significant error. The error from PB(n) on the 
other hand is only due to measurement noise. For this reason, it 
should turn out that the identification method presented v;ill be most 
successful when the boundary buses have relatively large power injec- 
tions (generators or loads) which supply m.uch of the known variation. 

Or in other words, it would be nice to have good, strong input signals. 
In actuality IS may be of relatively large size so that contributions 
of ^(n) to PEQl (n) = API Pl(n) from, far across the network may be of 
questionable value and accuracy due to the linear approximation of AQI . 
Since intuitively, buses in IS closest to BS will have the most affect 
on PEQl(n) and involve the least error, it m.ay be practical to use only 
part of the contribution of Pl(n) to PEQl (n) and lump the rest with 
PSQX (n). 



7.2 kodel Verification 

Although contrary argum.ents could be made, the boundary bus 
impedance model is probably the better of the tvro m^odels since it does 
not in\rolve estiir.ating a significant bias. However, the tie line power 



flow irodel is still valid for verification purposes. For the boundary 
bus irr-pedance model from equation (4.3) 



N 



5=^2 c(n)y(n)y' (n) 



n=l 



N 



= i 2 c(n)ZM ^X (n) ^QX ' (n)ZBB 



/X 



N 



(7.1) 



n=l 



A /\ 

If R is prerultiplied and post multiplied by ZB3 then an estim.ate of 



A /N, -1 

Bp£QX = ZB3 ZB3 ^ 



the covariance of PEOX(n) is 

(7.2) 

I'bw if measurements of PTL (n) were taken at the same times as measure- 
m:ents of D3(n) and FFQl (n) + ^(n), and if the estimate of YEQX is 

YEQX = ^ - YEQI - Y33 

/N A 

then equation (7.3) evaluated for YEQX should be close to SpgQx of 
equation (7.2). 

N 



\ 2 c(n) FEQX (n) PEOX ' (n) 



n=l 



I 2 -r[£ll(o) - (YTL + YEQX )DB(n)'l 
n=l 

X [fTI'(o) - (YTL + YEQX )D3(n)V c(n)"^ 



(7.3) 



If the two were not close, then the tie line power flow model would be 
indicating a lack of validity for the estimate of ZBB from the boundary 
bus impedance model. The tie line power flow model could also be used 
"on-line" once YEQ X had been identified to provide a continuing check cf 
the validity of the identified equivalent system. Should the m.easured 
value of PTL(n) vary from ( YTL + YEQX )D3 (n) an amount not commensurate 
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I 
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with y -PEQX’ vould indicate a change in XS had taken place. 

The change of course could either be in transmission line status 
affecting YEQX or a significant change in ^(n) affecting KQX(n). 



7.3 EstimatiPj? Equivalent Power Injections 

Once YE OX has been identified, the actual values of P EQX (t^^^ could 
be estimated. With enough m.easureir.ents it might be possible to develop 
an approximate daily pattern for PEQX (t^ so that they could be used as 
pseudo measurements [83 for the purpose of predicting a future power 
flow situation. If such were the case, the equivalent system could be 
extended in use from^ predicting changes in CS during which the change 
FEQX(n) would be small to use in predicting actual line power flow's for 
given power injections in CS. 



7.x Inputs V'hich Decrease Identification Error 

It has been assumed that the identification would be prim.arily 
passive in nature (i.e. by "listening" to the system) so that any power 
changes or voltage angle changes used to identify the equivalent system 
would be due to changes in the system's load demands. It was also 
argued that 5 ® reasonable approximation to the noise 

covariance where c(n) is inversely proportional to the time, t^^ - t^, 
between measurements . However, any large, relatively fast changes in 
OS v;ould provide an input, output measurement set with relatively little 
noise (FEOX(n)). These changes could be planned or accidental. If 



planned changes were used, it would be desirable to obtain a m.aximum of 
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invorir.ation with a miniinxim of change. For this case it might be 
possible to formulate the problem in terms of minimizing a cost 
function subject to a constraint such as finding the inputs 
u(n) (n=l,2,. . . which minimizes the cost 



7 .5 Use of the Equivalent System ^Jhen the F-xternal System is Known 

Although the linear load flow miodel was used to solve the problem 
of finding equivalent systems with no prior knov/ledge of the external 
systems, there is still an advantage to using the model when enough is 
known about the external system so that, for instance YEQX and PEQX 
could be found analytically for a particular system. The linear model 
clears the smoke so to speak, so that a good approximate picture can 
be easily obtained. For example AQX (and AQI ) gives a general overview 
of how power is routed by the network while a comparison of YEQX with 
YB3 shows at a glance the relative coupling between the two networks. 




subject to 



z(n) = H u(n) + y(n) (n = 1,2,...,N) 



where f^ and f^ are scalar functions 
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7 .6 Identification of the Nonlinear Kodel 

Use of the linear model may be sufficient to predict line power 
flows and voltage angles. However, in order to predict the effect of 
XS on voltage magnitudes within OS, the full nonlinear model is neces- 
sary, developed in Appendix C are the equations relating complex bus 
voltages to complex bus power and current injections. The computation 
required and complexity of identifying the nonlinear equivalent model 
would likely be very much mere than for the linear model. The increased 
accuracy in predicting real power flows would most likely be minimal, 
but the linear model cannot predict voltage changes due to the struc- 



ture of XS. 
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8. COrCLUSIOIS 
^.1 Areas for Further Study 

The tie line power flow model produces a biased estimate of the 
equivalent model parameters, however, if an unbiased estimate could be 
obtained so that accuracy was comparable to the boundary bus impedance 
model, there would probably be an overwhelming advantage to using the 
tie line power flow model. One of the most important reasons is that 
it uses variables from a local region only. PTL (t) is directly 
measurable and in order to find OT(t) state estimation would be required 
only in a small area around Boundary System^ (Power coming in from the 
rest of Own System could be treated the same as PTL (t^ for purposes of 
stats estimation.) Also the linear tie line power flow m^odel appears 
to be less of an approximation to its nonlinear analog (see Appendix C) 
than the linear boundary bus impedance model is to its nonlinear analog. 
Therefore, it is suggested that a study to find a method for eliminating 
bias in the tie line power flow model would be worthwhile. 

The simulation carried out used a linear load flow model for the 
puroose of generating data used to identify the equivalent model. This 
simplified the simulation and it also served to eliminate errors which 
would be caused by variations from ideal conditions (lew line resistance 
and voltage magnitude nearly constant). A more realistic simulation 
should be made to find how data from a nonlinear model affects accuracy 
of the identified linear model parameters. Also the net power flow summed 
over all the tie lines is usually held to a scheduled value. The effect 
of this on the identification should also be studied. 

As mentioned previously, the identification accuracy could be 
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increased by varying the generation pattern for the purpose of probing 
the external system. However it would be necessary to maximize inform.a- 
tion and minimize the variation required. It is felt that one step in 
the proper direction m.ight be to move one generator at a time since this 
would eliminate any interference among the inputs. A look at the boun- 
dary bus impedance model written in term.s of the general m.odel 
z(n) = H u(n) + v(n) shows that if elements of u(n) were all zero except 
for element i, then the i^^ column of K would be ^^^zCn) - v(n)]. 

A study into this problem should result in a more sophisticated solution 
to a rather com.plex problem. 

An obvious extension of identification of the linear model is the 
identification of the full nonlinear miodel. It is felt that the com- 
plexities involved would be an order of magnitude above that required 
for identification of the linear miodel. The accuracy gained in predicting 
voltage phase angle changes would probably be sm.all, but voltage m;agnitude 
changes which have been ignored, would be available from a nonlinear 
m.odel. If a search type of solution were used, then of course the 
solution to the linear model identification would give a starting point. 



H.2 Surmary 

The problem, of predicting how an electric power network will interact 
with other networks has been solved by reducing the entire outside world 
down to a relatively sm.all equivalent m.odel. The problem, is com.plicated 
by the facts that (1.) no knowledge of parameters or variables of the 
outside world can be assumed, (2.) identification m.ust be prim.arily 
passive, and (3.) "signal to noise ratios'' can be expected to be about 1:1. 
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A linear r^odel for the overall system vas developed and its use 
justified with respect to the more conventional nonlinear model by con- 
sidering accuracy versus computation requirements, insight, and com- 
plexity of the identification approach required. Two models were then 
proposed with each having certain advantages and shortcomings. The tie 
line power flow model uses data which is relatively accurate and locally 
available, but it produces a biased estimate. On the other hand the 
boundary bus impedance model produces an unbiased estim.ate but it 
requires data from all over the system and also requires manipulation 
of the equations for own system. It was argued that changes which occvr 
in power and bus voltage angle between two times should be used with the 
hope of obtaining a zero mean disturbance and also to avoid the }2se of 
matrices whose elements have magnitudes very large com.pared to the 
determinant of the matrix. Since both models were in the general form 
of a linear, static, input-output relation with additive disturbances, 
a solution to the general problem was found under assumed conditions. 

Then the solution was applied to the two models and it was shown how 
the tie line power flow m.odel produces a biased estimate. The simulation 
made to test the r.ethod was discussed shewing the excellent correspon- 
dence between theory and simulation results. Finally, practical aspects 
of the problem were briefly discussed including model verification, 
estimation of the disturbances, system probing, and use of the linear 
model reduction m:ethod for purposes other than identification. 

This method is proposed as a m^eans of modeling the outside world 
and finding the parameters of that model.. The actual solution is in a 
very simple form and could probably be implemented as an extension of 
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state estimation techniques. It is by no means the ideal answer to how 
the outside world should be modeled, but it does allow a prediction of 



how the outside world will interact with the system of interest 
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APPEIDK A. N:)M.IIEAR U^AD FLW EQUATIOJS 



For the pi transmission lire model shown in Figiire ?.l, the current 
flowing into the line from bus i is 



>'=iic ^ik 



(A.l) 



when the line is connected between buses i and k. is the shunt 

capacitive admittance, is the series admittance, and and Ej^ are 
the voltages at buses i and k respectively. All are complex quantities 
and 



— . . *^ik 

y^k = jy^ik = T 



-j0 



^ik = ^ik ® 



ik 



R., + jX., 
ik ik 



E. = E . e 
i 1 






^k ' ^k ® 



jO 



k 



(A. 2) 

(A.3) 

(A.4) 
(A. 5) 



where ^ik’ ^^ik resistance, reactance, and capacitive 

susceptance of the transmission line between buses i and k. and 

are the voltage phase angles at buses i and k as measured with respect 
to a reference. Equation (A.l) can be rewritten 



ttjj, = E.(y= 



ik 



^ik^ " ^ik^k 



(A. 6) 



By Kirchhoff's current law, the sum of currents entering the tran^m.ission 
lines connected to bus i must equal the current being injected into the 
bus by generators or loads (positive for generators and negative for 
loads). is the current injected into bus i. 
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? ^ik ^ik^ “ ^ ^ik^k 



(A. 7) 



Let a complex bus admittance matrix whose diagonal element in 

row i is the sum of shunt and series admittances of all lines connected 
to bus i 

-oe,, 



^-bushi ■ ^ik^ ^ii ® 



11 



(A. 8) 



and whose off diagonal ik and ki elements are the negative value of the 
series admittance of the line between buses i ard k (or zero if there is 
no line between buses i and k). 



^-bushk "^ik ^ik ~ ^ik ® 



"'’®ik 



(A. 9) 



-bus ® syTiiTietric matrix. Then for com.plex vectors -bus 

whose i^^ elements are the current injected into bus i and the voltage 
at bus i respectively, equation (A. 7) can be written in matrix form as 



■^us “bus “bus 



(A. 10) 



When Kirchhoff’s current law is applied to the entire network, it is 

obvious that for a system with K+1 buses, the current injections into 

bh 

K of the buses determine the current injected at the (K+l)^^ bus. 
Similarly, the voltage angles are relative to one another so one bus 
must be designated the reference bus for voltage angle whereas the 
reference for voltage m.agnitude is ground. 

Real and reactive power injected into bus i, + jQ^, in term^ of 



the current injected and bus voltage is 

P, + jQ, = E. V 
Pi - oQi - S/ii 



(A. 11) 



or 



(A. 12) 



where is the complex conjugate of voltage at bus 1. Using 
equation (A. 7) along with (A. 8) and (A. 9) where is defined, 
equation (A. 19) is 



P. 

1 




K+1 




ik 



> g ® ^ik ® 



-j^ik 



P. 

1 



k=l 



(A. 13) 



(A. 14) 



The load flow equation, (A. 14), can be used for Tietworks consisting of 
transmission lines and fixed tap transformers. However for other cir- 
cumstances such as voltage controlled buses, tap changing transformers, 
and phase shifting tran^ f orm.ers the equation must be modified. Stagg 
and El-Abiad CIO] is one reference for these arid other variations. 
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APFEJDK 3. ITiAD FLaJ EQUATI0I5 

The linear approximation for real power flow in a transmission 
line between bus i and bus k is 



tance magnitude. This approximation is subject to the conditions 

1. D. - D, is small. 

1 k 

2. The bus voltage magnitudes do not vary much from their 
nominal values. 

3. The ratio of transmission line resistance to reactance is 
small. 

For convenience it will be assumed that = 1.0 (or equivalently 

y^^ is normalized for E^ and E^^) so that equation (B.l) is 



Assvaning no losses within the bus, by continuity of power the power 



loads) is equal to the algebraic sum of the power flows into the trans- 




(3.1) 



where and are the voltage phase angles at buses i and k as 
measured with respect to a reference, and are the nominal values 
of the bus voltage magnitudes, and is the transmission line admit- 




(B.2) 



being injected into a bus (positive for generators and negative for 



mission lines corrected to that bus. Let be the power injection at 
bus i, then 



P. 

1 




(B.3) 



Define P to be a vector whose i 



th 



element is the power injected into bus 



i, and define D to be a vector of corresponding bus voltage angles. 



I 



^ip 
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Then relation (B. 3 ) becomes 
P = y D 



(3.4) 



E;/ comparison with (B. 3 ) it can be seen that the elements of Y are such 

that the diagonal eler.ent of row i is the sum of admittance magnitudes 

of the transmission lines which are connected to bus i, Y.. = 2 y.i • 

k 

The off diagonal elements ik and ki are negative value of the admittance 
m.agnitude of the line between bus i and bus k, Y^j^ = Y^^^ = 

Therefore the miatrix Y is of the form 



Y = 



2 

k 

-y 

-y 



yik 


"^12 


!>> 

1 


12 




-Y 23 


13 


-Y 23 





(3.5) 



Because each bus is usually only connected to a few of all the other 
buses, Y is normally a sparce matrix. It is readily observable that 
the sum of all elements of any row or column of Y is zero. The bus 
voltage phase angles must be measured with respect to a common reference 
so one bus must be designated as the reference bus where voltage angle 
is specified. Also, because the linear model is lossless, by continuity 
of power, the algebraic sum of elements of P is zero. Hence one row of 
equation (3.4) is redundant. If bus i is chosen as reference for vol- 
tage phase angle, then let vector D with elem.ent i deleted, 

let P, be the vector P with element i deleted, and let Y. be the 
-bus - -bus 

matrix Y with row i and column i deleted. Then equation (B.4) is 

P, = y, D, 

-bus -bus -bus 



(B. 6 ) 



I 
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P, and D, can be separated into vectors PI and DI for Internal 
“bus -bus ‘ — — 

System (IS), K and ^ for Boundary System (BS), and K and DX for 
External System (XS). Similarly separating into submatrices 



PI 




YII 


Y13 


IE. 






FB 


= 


YB» 


YB3 


YBX 




D3 






YIX» 


YBX ’ 


YXX 




BX 


^ __ 








— 




— 



(3.7) 



However by definition there are no transmission lines between buses of 
]S and XS so XJC = 0. 



II 




YII 


YB 


0 








FB 


= 


YB» 


YBB 


YBX 




D3 


(B.8) 


PX 




0 


YBX* 


YXX 




DX 




— — 












- 





The inverse of the bus admittance matrix is the bus im.pedance 



matrix Z, such that 
-bus 



D. 



-bus 



^us -1 



•bus 



(B.9) 



or 


HI 






D3 


= 




DX 











ZII 

ZI3 * 

ZK' 



ZIB 

ZB3 

ZBX» 



ZIX 



ZBX 

ZXX 



II 

PE 

PX 



Z, in term:s of submatrices of Y, can be found bv solving 
-bus -bus 



^us 


Y, 

-ous 


= I for 


-bus 


wh 


ere J is the 


identity inatrix. 






















— 




ZII 


ZB 


ZIX 




YII YB 


0 




OI 


0 




ZB* 


ZBB 


ZBX 




YB* YBB 


YBX 




0 I 


0 




ZIX* 


ZBX* 


m 




0 YBX* 


YXX 




lo 

lo 


I 












— 











(B . 10) 



(B.ll) 
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Kultiplying the middle row of Z, bv the left column of Y. (here 
^ ^ -^us " -ous 

rows and columns will refer to row and column submatrices.) 

ZB' YT^+ ^ YB' = 0 (B.12) 

and solving for ZB ' 

ZB ' = - Z3B YB ' YII~^ (B.I3) 



Kultiplying the middle row of by the right column of 

solving for Z3X 

ZB3 YBX + m Y^ = 0 

Z3X = - Z33 YBX YXX~^ 

Multiplying the middle row of Z^^^ by the middle column of 
solving for Z33 with the use of (3.I3) and (3.15) 

ZB ' YB + Z^ Y^ + ZBX YBX' = I 

- ZB3 YB ' Yir^ YB + Z3B - ZB3 YM YXX'^YBX » = I 

ZB3 = (-YB' Yir^ YB + XB3 - XiX XXX"^ YBX')~^ 



(B.I4) 

(B.I5) 



(B.16) 

(B.17) 

(B.18) 



Multiplying the top row of Z^^^ by the left column of X^^g and solving 
for ZII with the use of (3.13) 



ZII YII + ZB YB' = I 



ZII = YII“^ - ZB YB' YII"^ 



ZII = YII"^ + YII"^ YB ZEB YB' YIl"^ 



(B.I9) 

(3.20) 

(3.21) 



Multiplying the top row of Z^^^ by the right column of X-j^^g and solviTig 
for ZIX with the use of (B.I3) again 

ZB YBX + ^ XM = 0 (3.22) 



i 



i 

I 



L 



9 



I 
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ZIX = - ^ YXX ~ (B.?3) 

ZK = Yir H^IB ZH3 Y3X YXX'^ (3.24) 

Let matrices YSOI . YEQX. AQI . and ^ be defined by 

YSQI = -YB» YII"^ YB (3.25) 

YEQX = - YBX YBX * (3.26) 

AQI = - YIB ' YII"^ (3.27) 

AQX = - Y3X Y:« ~^ (3.28) 

Then equations for ZE3 . ZII . ZI3 . ZIX . and Z3X can be written 

Z3B = ( YEQI + ^ + YEQIC )"^ (3.29) 

ZII = YII"^ + AOI' ZBB Agl (3.30) 

ZB = A2I' ZBB (3.31) 

ZIX = AOI» ^ AQX (3.32) 

Z3X = ZBB A'QX (3.33) 

Using (B.IC), the equations for DI and DB are 







ZII 


ZB~ 




II 




2IX 




= 










+ 




D3 




ZB' 










Z3X 



or using (3.29) through (B.30) 
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+ 




ZHP AQX PX 



( 3 . 35 ) 



From equation 



(3.35) the key role of ZB3 is obvious. 
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AFPEjDIX C. rD^JlI?E;.R MODEL IDE M I? IC AT ION 



To avoid notation clutter in this appendix, all voltage and current 
vectors and admittance matrices will be understood to be complex vectors 
and matrices even though there is no bar over the symbol (for a complex 
quantity) or bar under the sjTnbol (for a vector or matrix). For example 
E and I are complex vectors. 

The complex matrix equation for relating bus voltages to bus cur- 
rent injections is 

I = Y, E, (C.l 

DUS bus bus 



I^^^ and can be divided into separate vectors for Internal System 

(ll and El respectively), for Boundary System; (IB and E3), and for 
External System (H and EX). Dividing into the corresponding 

subm.atrices and noting that by definition of IS, BS , and XS there is no 
immediate coupling between IS and XS, equation (C.l) can be written 



* II 




YII YIB 0 




El 


IB 


= 


YB' YBB YBX 




EB 


K 




0 YBX' YXX 




EX 


_ 




- 




— 



The only coupling betv/een CS and XS is due to YBB which can be separated 



into two matrices. That belonging to OS will be designated YBOS , and 
that belonging to XS will be designated YTL, so that YEB = YBOS + YTL. 
separating equation (C.2) 
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II 

IB 

K 



YII yiB 0 

YIB' YBCS 0 

0 0 0 



El 

EB 

EX 



0 

0 

0 



0 0 
YTL YBX 

YBX ' YXX 



El 

EB 

EX 



(C.3) 



Removing the equation for OS 



II 






YII YIB 



YB YBOS 



El 



E3 



YTL EB + YBX EX 



(C.4) 



The quantity YTL EB + YBX EX is a vector of line current leaving the 
buses for BS over the tie lines between BS and XS. Designating this 
the tie line vector TIL 



II 




E-TTL 





YII 


YB 


YE' 


YBOS 



EB 

TTL = YTL EB + YBX EX 

Removing the bottom row from equation (0,3) 
IX = YBX' EB + YXX EX 
and solving for EX 



EX = -Y:(X"^ YBX' EB + YXX"^ IX 



Substituting equation (C.8) into (0.6) 

TTL = YTL EB - YBX YXX"^ YBX' EB + ITX YXX"^ K 



(0.5) 



( 0 . 6 ) 



(0.7) 



( 0 . 8 ) 



(0.9) 



Define an equivalent external bus admittance matrix YEQX and an equi- 
valent external complex bus current injection vector lEQX by 
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YEQX = -YBX YXX"^ Y3X» 


(C.IO) 




EQX = AQX IX 


(C.ll) 


where 


AQX = -YBX YXX"^ 


(C.12) 



Then equation (C.6) becones 

TTL = (YTL + YEQX)EB - lEQX (C.13) 



and if equation (C.I 3 ) is substituted into equation (C. 5 ), since 
YBB = YBOS + YTL 



II 




YII 


YB 




El 


B3+EQI 




_YB' 


(YBB+YEQX) 




EB 



In the sair.e vay, an equivalent complex bus admittance matrix and an 
equivalent complex bus current injection vector for E as seen by BS 
can be found. The top row of equation (C.I 4 ) is 

II = YII El + YIB EB (C.I 5 ) 

Solving for El 

El = -YIl"^ YIB EB + YII"^ II (C.I 6 ) 

and substituting into the bottom row of equation (C.14) 

B + EQI = -YB' Yir^ YB EB + YB» YII"^ II + (YBB + 1EQX)EB (C.17) 



Defining an equivalent internal complex bus admittance matrix YEQI, and 
an equivalent internal complex bus current injection vector lEQI by 





YEQI = -YB» YII"^ YB 


(C.18) 




BQI = AQI II 


(C.19) 


where 


AQI = -YB' YII"^ 


(C.20) 
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equation (C.17) becomes 

lEQI + IB + lEQX = (YEQI + YBB + YEQX)EB 



(C.21) 



Just as YBB was divided into YTL and Y3CS, so also Y30S can be divided 
into that belonging to IS designated as YBES, with Y3CS = YBIS + YBBS, 
Then equation (C,2l) is 

EQI + IE + EQX = (YBIS + YEQl)EB + YBBS E3 + (YTL + YEQX)EB (C.22) 



Thus far equations (C.5), (C.I3), (C.I 4 ), and (C.2l) are the 
pertinent equations in terms of bus current injections. To place them 
in terms of bus power injections let a diagonal matrix whose elements 
are the complex conjugate values of bus voltages be 




E* 0 0 — 

0 E* 0 — 

0 0 eJ — 



(I'bte that [El is a diagonal m.atrix whose diagonal elements are iden- 
tical to the elem.ents of the vector E.) Also, let PI and QI, PB and 
Q3 , and PX and QX be the real and reactive bus pov;er injections of 
IS, 3S, and XS. Then because for any bus i, ~ 



PI - jQI = [eI*1 II 


(C. 23 ) 


PB - jQB = [EB"] IB 


(C.24) 


PX - jQX = [eX*J EC 


(C.25) 


FTL - jQTL = [eb'*] TTL 


(C.25a) 



or 



! 



• f 



I 
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II = [ei''] (pi - jQl) (C.26) 

DC = [ex"] (PX - jQX) (C.27) 

Multiplying equation (C.2l) by [EB ] 

[E3*](IEQI + IB + lEQX) = [eB*](YEQI + YEB + Y£n;C)EB (C.28) 

Define equivalent internal and external bus power injection vectors by 
PEQI - jQEQI = [eB*] IEQI = [eb'""] AQl[Ef] (PI - jQl) (C.2S) 

PEQX - jQEqX = [eB^] IEQX = [eb‘] AQX[EX*] (PX - jQX) (C. 30 ) 



Then (C.2S) becomes 

(PEQI - jQEQI) + (FB - jQB) + (PEQX - jQEQX) = [eb”] (YEQI + YBB + YEQX)EB 

(C.3I) 



• Sl T-rARY 

For the following a bar over a symbol denotes a complex quantity 
v/hile matrices and vectors are denoted by underlined symbols. Using 
equations (C. 5 ), (C.23)> (C.24), and (C,25a) the nonlinear load flow 
equation for CS is 



PI - m 




[ei"] 


0 




YII 


YB 






(PB - OQB) - (PTE - -iQTL) 




0 


[ra"] 




yb' 






H 



(C.32) 



which, if YEQX and PEQX - j QEQX are known, becomes 









I 



I 






Si 
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^PI - 


m 




\ 

1 

Ml 

yx 

C 1 


0 




YII YJ3 




El 


(ra - 


.10.3) + (PEOX - 1QEQX) 




0 






Yb' (Y3B+YEQX) 







(C.33) 



The nonlinear analog of the linear load flow tie line power flow model 
for identifying YEQX is 

PTL - jOTL = ( YIT. + YEQ X)S3 - ( FF.OX - i QSQX ) (C.34) 

and the nonlinear analog of the linear boundary bus impedance model is 



^( FSQI - .1 QSQI ) + (ra - + ( FEQX - j OEQX ^ 

= + Si + il 

Repeating the equations for current injections 



n 




YII 


YB 






S - SI 




IS' 


Si 




S 


s 


' 


Yii 


YB 




- 


B + BQX 




YB' 


(YBB+YEOX) 



ITL = (YTL + XE0X)EB - lEQX 
( lEQI + S) + = ( YEQI + ill + 



(C.35) 



(C.36) 



(C.3V) 



(C.3B) 



(C.39) 
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APPE ^D n D . K/.X KU1-' LIKE LIHCOD IDE !vT IF IC AT 10 N EQUAT 10 fS 



Given a sequence of input rneasurerrents u(l), u(?),..., u(N)> 
a sequence of correspon(5ing output rneasureir.ents z(l), z(2) , . . . , z(N), 
a set of confidence coefficients c(l), c(2),..., c(N) which indicate 
the relative confidence associated with each input/output set, and 
given the input /out put relation 



z(n) = H u(n) + v(n) 



(D.l) 



where v(n) is a noise vector having a Gaussian probability distribution 

1 



with zero mean and covariance 






R, it is desired to find the ma^ciuTUin 



likelihood estimates of K and R subject to the following conditions: 

1. The elements of H are independent of one another. 

(e.g. in general H / H') 

2. E ^v(n)^ =0 

3. E ^y(n)y' (rn)"^ = 0 for n / m 

4. E^y(n)y'(n)^ = - 

5. E [y(n) u' (n)^ = 0 

The probability distribution function for the Kxl vector y(n) is 



pfv(n)] 



r K K t" 2 -p Z' (ji)c(n)H"^v(n) 

= [{2iTf c‘^(n)lRl] ^ e ^ 



(D.2) 



The likelihood function for z(n) and u(n) given H and R where z(n) and 
u(n) are Kxl vectors, and H and R are KxK matrices is 



p[z(n),u(n); H,R] = [(2 tt)^ c"^(n)lRl] ^ gjrz(n) ,u(n) :H,R] 



where 



m 






I 



t 



9 
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J[z(n) ,u(n) :H,R] = |[z(n) - H u(n)] c (n)R"^[z(n) - H u(n)1 (D.4) 



Using to rcean the trace of matrix A and iteing the identity 

a’a = tr|^ a'"^ where a is a vector equation (D.4) can also be wrritten 

J[z(n) ,u(n) :H,R] = | tr |^c (n)R~^|^z (n) - H u(n)][z(n) - H uCn)")^ (D.5) 

let the sets Z, U, and V be such that Z = £z(l), z(2),..., z(N)^ > 

U = ^u(l), u(2),... , u(N)^, and Y = (l) , y (2) , . . . , y ( N)^ . 

Because y(n) is assumed to be uncorrelated in time (U^y (n)y ' (m)^ = 0 
for n / ra) and because y(n) is Gaussian, the joint probability distri- 
bution function for V is the product of the distribution functions for 
y(l), y(2),..., y(N). 



p(Y) = P[v(l)l p[v(2)1 . • • p[v(N)] 



N 



n=l 



1 -'5 S Y' (j^)c(n)R~^y(n) 



2 n=l 

e 



(D.6) 



The likelihood function of Z and U given K and R is then 



N 



P(Z,U:H,B) = ]X 

n=l 




-J(Z,U:H,R) 

e 



(D.7) 



where 



J(Z,U:H,R) 



I tr-^c(n)R ^[z(n) - H u(n)3[z(n) - H u(n)] ^(D.S) 



It is desired to find the values of H and R which maximize the likelihood 
function p(Z,U:H,R). These values of H and R are the maximm likelihood 

A A 

estimates H and R, If the logarithm of the likelihood function is m.axi- 



mized, then so will the function. The log likelihood function is 
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/n[p(Z,U:H,R)] = - | /h(2rr) f | g K /n c(n) 

n=l 

- I N/n[Rl - J(Z,U:H,R) (D.9) 



A new function can be defined 

N 

f(Z,U:H,R) = -2/n[p(Z,II:H,R)] - M /n(2rr) f ^ ^ c(n) 

n=l 



= N ^nfRl + 2J(Z,U:H,R) (D.IO) 



f(Z,U:H,R) only has terms which involve H or R, hence, minimizing 
f(Z,U:H,R) is equivalent to m.aximizing the likelihood function. 

To perform the necessary mathematical manipulations certain 
matrix gradient identities are used from Athans and Schweppe D-l. 

The matrix gradient for a scalar function of a KxK matrix X is defined 
for use here as 
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Two properties 





" 9f(X) 


^f(X) 




pxii 


PXki 


Cf(x)] = 


c9f(X) 


^»f(X) 






_ 


of trace of 


a matrix are 




tr{A} = 


tr ^A»^ 





tr 





(D.ll) 



(D.12) 



The identities used are 



= A 



fx 



(D.14) 












• I 



u 












I 



9 



I 






I 







I 
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p tr X B X'^ = 3 X' A + B' X' A' 



^ tr ^A X"^ = -X~^ B A X‘^ 



/n|X( = X'^ 



(D.15) 

(D.16) 

(D.17) 



When i'(Z,U,H,R) is at its minimum value, the matrix gradients of 
f(Z,U:H,R) with respect to H and R v;ill be zero. Then the maxirr;Ui:i like- 

A A 

lihood estimates H and R will be the solutions for K and R of 



— f(Z,U:H,R) = 0 



and 



for 



and 






c?R 



f(Z,U:H,R) = 0 



f(Z,U;H,R) = 2T/n|R| + 2J(Z,U:H,R) 



(D.18) 

(D.19) 

(D.20) 



N 



J(Z,B:H,R) = 5 Z tr ^c(n)R*^[z(n) -Hu(n)][z(n) 

n=l ‘■ 

First, the gradient with respect to H is 



H u(n)] '|(D. 



21 ) 



c2H 



f(Z,U;H,R) = “ 2J(Z,U;H,R) = 0 



c?K 



^ ^ tr ^c(n)R~\z(n) - H u(n)l[z(n) - H u(n)l^ 
- n=l ^ 

tr^c(n)R ^z(n)z^(n) - c(n)R ^z(n)u'(n)K' 
- c(n)R~^H u('n)z' (n) + c(n)R u(n)u' (n)H'J. 



<P 

n=l - 



(D.22) 



I 



I 



1 
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Using equations (D.12) and (D.I 3 ) 



tr (n)R~^z(n)u’ (n)K'^ = tr |"c (n)H u(n)z ' (n)R“^J’ 

= tr (n)u(n)z' (n)R 

tr (n)R ^ u(n)z ' (n)^ = tr -^c (n)u(n)z ’ (n)R 



Then equation (B.22) is 
N 



P = ^ tr ^-2c (n)u(n)z ’ (n)R ^ + c(n)R u(n)u* (n)H'^ 



and using equations (D.I 4 ) and (D.I 5 ) 



N 



0=2 I“2c (n)u(n)z ' (n)R ^ + 2c (n)u (n)u' (n)H'R 
n=l 

r N N ^ _ 

= 2 I- 2 c(n)u(n)z’ (n) + 2 c (n)u(n)u' (n)H’l R""^ 
*■ n=l n=l ^ 



,-l 



Then unless R =0 (which would mean infinite noise covariance) 



N K 

2 c (n)u(n)z' (n) = 2 c (n)u(n)u' (n)H’ 
n=l n=l 

N N 

K 2 c(n)n(n)u' (n) = 2 c (n)z (n)u’ (n) 
n=l n=l 



and the value of H v/hich maximizes p(2,U:H,R) is 



N T 


r N 1 


2 c(n)z(n)u» (n)l 


1 2 c(n)u(n)u»(n) 


,n=l -1 


ln=l 



Taking the matrix gradient of f(Z,U;H,R) with respect to R 



(D.23) 

(D.24) 

(D.25) 

(D. 26 ) 

(D.27) 

(D.28) 

(D.29) 
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«?R 



f(Z,U:H,R) = 0 



= N 



/nlRl + ^ tr |*c(n)R"^[z(n) - H u(n)] [z(n) - H u(n)]^ 

- n=l ^ 



(D.30) 



Using equations (D.16) and (D.17) 



N 



0 = iR ^ c(n)R ^l^z(n) - H u(n)]^z(n) - H u(n)^ R ^ 



n=l 



= m - y, 2 c(n)[z(n) - H u(n)][z(n) - H u(n)]'^R’^ 



(D.3I) 



Again unless R" = 0, using equation (B.29) the value of R which mini- 



mizes p(Z,U:K,R) is 






R = 4 2 c(n)[z(n) - H u(n)][z(n) - H u(n)] 



(D.32) 



n=l 



Hence H and R or equations (D.29) and (D.32) are the maxLmira likelihood 



solutions to the problem. 
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